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Particle-in-cell (PIC) is the most used algorithm to perform self-consistent tracking of intense 
charged particle beams. It is based on depositing macro-particles on a grid, and subsequently 
solving on it the Poisson equation. It is well known that PIC algorithms occupy intrinsic limitations 
as they introduce numerical noise. Although not significant for short-term tracking, this becomes 
important in simulations for circular machines over millions of turns as it may induce artificial 
diffusion of the beam. 

In this work, we present a modeling of numerical noise induced by PIC algorithms, and discuss its 
influence on particle dynamics. The combined effect of particle tracking and noise created by PIC 
algorithms leads to correlated or decorrelated numerical noise. For decorrelated numerical noise 
we derive a scaling law for the simulation parameters, allowing an estimate of artificial emittance 
growth. Lastly, the effect of correlated numerical noise is discussed, and a mitigation strategy is 
proposed. 

PACS numbers: 29.27.Bd, 29.27.-a, 41.75.-i, 52.59.Sa, 


I. INTRODUCTION 

It is well known that in operational scenarios requir¬ 
ing long-term storage, sources of noise in the machine 
can lead to detrimental effects on the beam. Exam¬ 
ples include RF voltage noise, beam-beam interaction, 
ground motion, and feedback noise, which all have been 
subject of studies to mitigate or exploit their impact, as 
reported in Refs. [1-6] 

Recently, a similar concern has raised for numeri¬ 
cal noise in self-consistent simulation of high intensity 
beams. In fact, the study of space charge effects via 
self-consistent simulations has become important with 
the advent of new projects, like the future SISlOO syn¬ 
chrotron of the FAIR project [7], and the LIU project [8] 
for the CERN accelerator complex. In these projects 
some scenarios require the storage of a high intensity 
bunched beam for seconds. 

The simulation of these operational scenarios signif¬ 
icantly raises the computational challenges, provided 
that they concern the prediction of beam loss, or emit¬ 
tance growth. In order to avoid any artificial noise, 
the first simulation studies were limited to frozen mod¬ 
els, which rely on noise-free tracking schemes [9, 10]. 
There, the space charge is computed assuming the beam 
remains frozen, therefore allowing an analytic descrip¬ 
tion of the space charge force. This approach surely 
avoids artificial emittance growth, but unfortunately at 
the expenses of self-consistency, which becomes relevant 
in certain scenarios e.g. large beam loss, or large beam 
core growth. 

On the other hand, PIC simulations of high intensity 
beams allow a self-consistent treatment, but require an 
understanding of the origins and propagation of its nu¬ 
merical noise. Hence, simulation parameters can be set 
adequately in order to mitigate noise induced artifacts, 
while keeping the computational load feasible. 

The evaluation of PIC induced artifacts on particle 
tracking has been studied by several authors from sev¬ 
eral points of view, e.g. in the context of a Eokker- 


Planck approach, as reported in Refs. [11-13]. 

Differently from previous approaches, we discuss the 
effect of PIC induced noise by following the integration 
method of a typical beam dynamics code. We start 
with an analysis of the noise due to the PIC algorithm, 
and discuss how it propagates via the beam dynamics 
integration. Eollowing this approach, the dependence 
of rms-emittance growth on simulation parameters is 
derived. It is also found that the particle’s rotation 
in phase space creates correlations in numerical noise, 
which enhances artificial emittance growth. 

Our results are of interest for simulations of millions 
of turns in circular machines, where the high intensity 
requires a self-consistent modeling, hence the use of PIC 
algorithms, while at the same time the control of artifi¬ 
cial emittance growth. 

This paper is organized as follows. In Section II we 
study the origins and properties of static PIC noise. 
Then, in Sec. Ill, we model the effect of random PIC 
noise on the dynamics of a single particle, while the 
effect on the whole ensemble of particles is treated in 
Sec. IV. In Section V we introduce the concept of cor¬ 
related and decorrelated numerical noise. We conclude 
our studies with a summary in Sec. VI. The effect 
of the particle’s rotation in phase space on numerical 
noise is discussed in Appendix A, whereas the conse¬ 
quences of the periodic random walk on the diffusion 
of macro-particles is treated in Appendix B. Finally, 
in Appendix C, we discuss the impact of the strength 
of space charge forces on the excitation of stochastic 
resonances. 


II. SCALING LAW FOR ELECTRIC FIELD 
FLUCTUATIONS 

Before treating the propagation of numerical noise, 
we study how the PIC scheme is causing numerical noise 
in a static scenario. We start by considering a Gaus¬ 
sian beam, i.e. the particle’s phase space positions are 
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defined randomly according to a Gaussian probability 
density function (p.d.f.), given by 


f{x,x',y,y') = 


' J-V _L 2/ 


47r^(J x^x' 


( 1 ) 


where cTxj'Jx',<T y,ayi are the transverse standard devi¬ 
ations. In the following discussion, we consider round 
Gaussian beams with := 2 ct^ = 2ay the standard 
deviation for the radial coordinate = x'^ + y'^. The 
associated radial space charge field at a longitudinal po¬ 
sition Sj is given by 


Er{x,y,Sj) 


Eg{Sj) 
\/x‘^ +y‘^ 



( 2 ) 


where Eg{sj) is proportional to the gradient of the elec¬ 
tric field at a: = y = 0. 

The transverse electric field of a coasting beam, that 
determines the self-consistent space charge forces, is 
computed in long-term simulations most efficiently us¬ 
ing a two-dimensional particle-in-cell (PIG) scheme. To 
achieve a high computational efficiency, the Poisson 
equation in the transverse plane is solved on a finite 
set of No X No grid points and the particle distribu¬ 
tion is approximated by Nm macro-particles. Since the 
number of macro-particles Nm is much smaller than the 
number of physical particles Np, an artificial granular¬ 
ity of the distribution is introduced, which is one source 
of numerical noise. 

To study the impact of the simulation parameters 
on numerical noise, we adopt the random start tech¬ 
nique (opposite to the quiet start). In the random start 
approach the positions of Nm macro-particles are ran¬ 
domly initialized according to a p.d.f, as given by Eq. I 
for Gaussian beams. Therefore, the resulting electric 
field is not only determined by the p.d.f. itself, but 
varies according to the random initialization of macro¬ 
particles with an amplitude 5Ex. This amplitude of 
fluctuations SEx is identified with the standard devia¬ 
tion of the electric field for multiple random initializa¬ 
tions of macro-particles. 

Previous studies, e.g. Refs. [14-18], have found that 
the parameters of relevance for a noise analysis of a 
beam with a fixed size are Nq and Nm, since they fix 
the number of particles per cell. To find the scaling 
of SEx on these parameters, we fix one of them and 
vary systematically the remaining one. We find that 
the simulation data satisfies the scaling 

SExix = 0,y = 0,s = So) (X , (3) 

\Nm 

where the dependence on the grid resolution, given by 
the function gi^No), is determined by the specific inte¬ 
gration method of Poisson’s equation. For the multi¬ 
particle tracking library MICROMAP [19], that incor¬ 
porates a spectral method, we find that g{Na) ~ ^Nq 
for well resolved beams, which is used throughout the 
article. In particular, we find that the level of numerical 
noise can be kept constant, if we choose our simulation 
parameters such that 


= const. (4) 

Nm 


A dependence of the electric field fluctuation on x and 
y is created for non-constant particle distributions, as 
the number of macro-particles per cell varies spatially. 
For round Gaussian beams this is given by the square 
root of the p.d.f., as we find 


SEx{x,y,s) 


SEx,o^ 


e ^{x,y,s) . 


( 5 ) 


Here, 


SEx,o = 5Ex{Q,Q,s) ( 6 ) 

y \/ng,o 

is a factor to quantify the standard deviation of the 
electric field. It is obtained by calculating the standard 
deviation in the center of the beam, (5F3,(0,0, s), for 
multiple random start initializations for Nm,o macro¬ 
particles and Ng,o x Ng,o grid points. We normalize 
(5Fa;(0,0, s) to Ng,o and Nm,o, such that SEx^ is inde¬ 
pendent of these simulation parameters. 

The factor ^{x,y,s) describes the effect due to the bi¬ 
linear interpolation of the electric field in between grid 
points, which causes a grid texture. As an example, the 
standard deviation for No = 16 and Nm = 2000 and a 
round Gaussian beam is given in Fig. 1. The statistics 
is done with 500 random start initializations. 
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FIG. 1: (Color) Standard deviation of the transverse electric 
field, see Eq. 5, for a round Gaussian beam. The exponential 
spatial dependence of the field fluctuations is superimposed 
by a grid texture caused by the mesh of the PIC solver. The 
simulation parameters are given in the text. 

The scaling law, Eq. 5, is valid as long as the beam is 
reasonably resolved, which is the case for No > 16 grid 
points within [—2ax,2ax] (and resp. in [—2ay,2ay]) 
of the Gaussian particle distribution. The effects of a 
low mesh resolution were recently studied in Ref. [13], 
where it is shown that due to grid heating the emittance 
growth is artificially enhanced. In the studies presented 
in this paper, we consider sufficiently resolved particle 
distributions, as this is the preferential case in space 
charge simulations. 

Using the scheme presented in this section, we can 
easily obtain scaling laws for other particle distribu- 
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tions. In the following, we study the Kapchinski- 
Vladimirski (K - V) particle distribution [20], i.e. a 
spatially constant distribution of particles in the trans¬ 
verse plane, causing linear space charge forces inside the 
beam. For this distribution we find that the electric field 
fluctuations dEx{x,y,s = sq) and 6Ey(x,y,s = sq) in¬ 
side the beam are constant, since the uncertainty of a 
PIC solver depends on the number of macro-particles 
per cell. Therefore, the electric field fluctuations inside 
a K - V beam are given by 


SE^{x,y,s) 


SE,,o^ 



^{x,y,s). 


(7) 


In Figure 2, we show the spatial dependence of the 
electric field fluctuation SE^ix^y, s = sq) for a K - 
V distribution of size X = 2ax and a rms equivalent 
Gaussian beam [21, 22] of rms size ax- For this study 
we took 1000 samples of electric fields (random start), 
while the beam was well resolved with a fine grained 
mesh of Nq x Nq = 512 x 512 grid points. 



X - position 

FIG. 2: (Color) Electric field fluctuations 5Ex{x,y) for a 
K - V distribution of size X = 2ax and a rms equivalent 
Gaussian beam. Inside the beam, i.e x € [—X,X], the field 
fluctuations for the K - V distributions are approximately 
constant. 

As seen in Fig. 2, the maximum electric field fluctu¬ 
ations for the Gaussian beam are by a factor larger 
than for a rms equivalent K - V beam. This is due 
to the fact, that in the center of the Gaussian beam 
the density is the double with respect to a rms equiva¬ 
lent K - V beam. The number of physical particles in 
the center is therefore by a factor 2 larger for Gaussian 
beams, which enhances the electric field fluctuations. 

The results obtained in this section apply for the 
numerical noise in coasting beam simulations. How¬ 
ever, the same method can also be used to find the de¬ 
pendence of the fluctuations SEx{x,y, s) for a bunched 
beam in the 2.5D scheme, see Ref. [23]. 


III. NOISE AFFECTING THE SINGLE 
PARTICLE DYNAMICS 

In this section, we study the effect of a fluctuating 
force field on the dynamics of a single particle by use 


of the random walk theory. Later, these findings can 
be applied on a large ensemble of particles, such that 
conclusions on artificial diffusion of a full beam can be 
drawn. 


A. Random walk model 

A particle i at the longitudinal position s is described 
by its transverse phase space coordinates Xi,x^,yi,y[. 
Its dynamics in a linear lattice is given by the Hill equa¬ 
tions 


x'l - kxXi = 0 , y'l - kyyi = 0 , ( 8 ) 

with kx, ky the restoring force terms. In the following, 
only the x-plane is discussed, since the y-plane can be 
treated equivalently. 

A particle i with a charge q and mass mg of a high 
intensity beam is exposed to a force created by the elec¬ 
tric field Ex{xi,yi, s) of the beam. The electric field 
is calculated in the rest frame of the beam, which is 
the standard approach in beam dynamics with space 
charge [24], and is used in this form throughout all dis¬ 
cussions in this paper. 

The effect of the electric self-field on particle dynam¬ 
ics is modeled via space charge kicks [25], whose effect 
is to change the particle coordinates x': 


^iisj) x[{sj)^ As 




(9) 


with c the speed of light, and /3 and 7 the relativis¬ 
tic factors. For simplicity, the longitudinal distance 
between two consecutive space charge kicks, given by 
As = Sj+i — Sj, is kept constant during the tracking. 

By knowing position and charge state of all particles, 
the resulting space charge field on a particle located 
at Xi,yi,Sj, given by Exfi(xi,yi, sj), is calculated as a 
superposition of all Goulomb fields. We refer to this as 
the exact solution. However, if approximations are used 
in order to reduce the computational cost of the space 
charge field calculation, imprecisions are generated. See 
e.g. Refs. [14-18] for a review. 

In the following, the effect of electric field fluctuations 
is studied in a simplified mathematical model, where the 
standard deviation is subtracted or added randomly to 
the exact solution, i.e. 


Ex {Xi ^ yi^ Sj') — 

Ex,o(^Xi ^yij s -t- Zij 5 Ex {xi ^ yij Sj^. (70) 

Here, Zij = 1 or Zij = — 1 is defined randomly with 
equal probability for any particle index i and any lon¬ 
gitudinal position index j. This model yields an av¬ 
erage of Exfi(xi,yi, Sj) and a standard deviation of 
SEx{xi,yi, Sj) for a large number of electric field cal¬ 
culations. 

The fluctuations Zij5Ex{xi, yi, Sj) are assumed to be 
statistically independent of each other, i.e. for every 
longitudinal position Sj, where a space charge kick is 
applied, the fluctuation Zij5Ex{xi,yi, Sj) does not de¬ 
pend on the fluctuation at any other longitudinal posi¬ 
tion Sj+fe = Sj + k ■ As for any fc G N. Throughout this 











4 


article, we refer to this kind of fluctuations as decorre- 
lated numerical noise. 

Equation 8, 9, and 10 are the model for the evolution 
of a single particle affected by numerical noise in a high 
intensity beam using the standard integration with the 
one kick approximation. We propose that the combined 
effect of numerical noise and transport can be modeled 
via an effective map: 


M 


x'r{Sj) + As 
_ ( X^{Sj+i) \ 


Xi{sj ) 

qlEx.o(xi,yi,Sj)+ZijSEx{xi,yi,Si)] 

moc2/32^3 


M 


0 

Zij Ax (xj, yi , Sj ) 


( 11 ) 


with Ax'{xi,yi, Sj) the strength of the random kick, 
that is dehned by 


Ax'{xi, yi, Sj) = As 


q5Ex{xi,yi,Sj) 


( 12 ) 


The left hand side of Eq. 11 is the usual one step in¬ 
tegration, where M is the transport matrix of a parti¬ 
cle between two consecutive space charge kicks obtained 
from Eq. 8. The coordinates Xi(sj+i), x'(sj-)-i) are those 
of the particle transported noise-free with space charge 
from Sj to Sj+i. The repeated application of Eq. 11 
mixes the noise ZijAx'{xi,yi,Sj) to both planes x and 
x'. Later, we model this effect by introducing an inte¬ 
grated effective noise kick. 

To model the effect of decorrelated numerical noise on 
Nm macro-particles at Nk space charge kick positions, 
we define two sets of independent variables Zij , Z[j with 
1 < i < Nm , 1 < j < Nk , in which Zij = — 1 or Zij = 1 
(resp. Z[j = — 1 or Zb = 1) defined randomly with 
equal probability for any index i, j. As suggested in 
Appendix A, we describe the integrated effect of nu¬ 
merical noise on the phase space coordinates by using 
an effective model: 


M 


0 

ZijAx'(xi,yr,Sj) 




( Zijf3x(^Sj) ^ Ax {xi,yi, Sj) 

I z[j J 7i 


(13) 


The additional random elements account for the inde¬ 
pendence of both planes x and x', while the normaliza¬ 
tion factor l/-\/2 stems directly from the mathematical 
model presented in Appendix A. The effect of Equa¬ 
tion 13 is to create a random walk [26]: At each longi¬ 
tudinal position Sj , where a space charge kick is applied 
on single particles, a random kick is added to both phase 
space coordinates. 

We test this ansatz in the following numerical exper¬ 
iment. A test-particle is initialized at x = 0 and x' = 0, 
while Nm = 1000 macro-particles are distributed 

according to a matched Gaussian probability density 
function. The beam is tracked for 1000 turns in a 
constant focusing channel of length L = 1 m with a 
tune set to Qx Qy 0.064, while the space charge 
induced tune shift is set to AQx — AQy ~ —0.005. 
The two-dimensional Poisson solver uses a mesh of 
Ng X Ng = 64 X 64 grid points. The integration 
applies one space charge kick per turn, such that 


As = L. The tracking is repeated n = 10000 times 
for the same physical conditions, but with different 
random initializations of the particle distribution. Due 
to numerical noise the test particle is found at different 
final positions. In Figure 3, we plot the distribution of 
final positions of the test particle after 1000 turns. The 
distribution is normalized to the number of particles 
found in the center of the beam, which is the maximum 
of the distribution. 
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FIG. 3: (Color) Distribution of final positions of a test par¬ 
ticle initialized at x = 0, x' = 0 after tracking in a constant 
focusing channel for n = 10000 different random initializa¬ 
tions of macro-particles. The distribution is normalized to 
the number of particles found in the center. 


The results presented in Fig. 3 show that numerical 
noise affects both phase space coordinates x and Pxx' 
in equal measure for a sufficiently large number of space 
charge kicks. 


B. Application: Single particle emittance 

In order to use the random walk model, we have 
to specify the spacial dependence of the fluctuations 
SEx{xi, yi, Sj). As a simplified approach we consider 

xf + yf 

SEx{xi,yi,Sj) = 6Ex{0,0) e '^’’r , (14) 

with SEx{0,0, Sj) = 6Ex{0,0) the electric field fluctua¬ 
tion in the center of the beam, which are, only in this 
first approach, assumed to be identical for any longitu¬ 
dinal position Sj. Equation 14 can be considered as a 
simplified model for the electric field fluctuations of a 
Gaussian beam, as derived previously, see Eq. 5. 

We now apply the random walk model, as given by 
Eq. 13, on a single particle i, that is initialized at the 
phase space position Xi^, x) q and tracked by using Eq. 8 
and 9. Following the standard theory of random walks, 
the most probable position of the particle after a series 
of p random kicks, is given by the position for a noiseless 
tracking. Figure 3 confirms this, as the maximum of the 
distribution of final positions is in the center. After p 
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random kicks, the variances of the particle’s phase space 
positions Xi^p, x' ^ scale as 

(xlp) - {x^,pf oc p5El{{), 0), 
{xf^p)-{x[^p?^p5El{Q,Q). 

Each particle can be associated with a single particle 
emittance 


^SP,Xi — 'yxXi “t“ ‘^CkxXiX^ PxXj^ , (1^) 

with Twiss parameters ax,Px-,"fx- Here, we discuss the 
effect of decorrelated numerical noise on the initial sin¬ 
gle particle emittance esp,xi,o in the context of the ex¬ 
pectation value of esp,xi,pi i-e. the single particle emit¬ 
tance after p random kicks. Assuming that numerical 
noise is the dominant cause of emittance growth, we 
write the expectation value of £sp,xi,p as a variation of 
the initial single particle emittance: 

{£SP,Xi,p) — £SP,Xi,0 + ^£SP,Xi,p- (17) 

Using Equations 15 and 16, we find the scaling 

^£SP,xi,p oc pSEl{0, 0). (18) 

Thus, in average, the single particle emittance grows 
over many realizations of space charge kicks due to 
decorrelated random fluctuations in the space charge 
field. 


IV. NOISE AFFECTING A DISTRIBUTION 
OF PARTICLES 

In this section, we apply the random walk model 
on each particle of the beam. We then find a scaling 

I 


on simulation parameters for noise induced diffusion. 
As an application, we predict the artificial emittance 
growth for a self-consistent space charge simulation for 
the SISlOO synchrotron at FAIR [7]. 


A. Effect of numerical noise on the rms-emittance 


The standard deviation of electric field fluctuations 
SEx{x,y, s) is given by Eq. 5 and respectively Eq. 7 
for any particle at position x, y, s. Since this result 
applies to any particle, it enables us to study the effect 
of decorrelated numerical noise on the full ensemble of 
particles, and in particular on the rms-emittance. 

In the following, we make use of the standard defini¬ 
tion of the rms-emittance in the x-plane 



4(s) = (x 2 )(x' 2 )-(xxr = 


Nm Nm 





(19) 


while the discussion for the y-plane is similar. The sym¬ 
bol (•) denotes the average of the particle distribution. 
Let ex{sj+i) be the beam emittance resulting from the 
application of a space charge kick at Sj and transported 
by As to Sj+i. By using the effective transport map 
Eq. 11, and with the noise modeling of Eq. 13, we pre¬ 
dict £a;(sj-|_i) from the beam distribution at Sj. In fact, 
from Eq. 19 it follows 




Z'-Ax'(xi(sj),yi(sj) 


Ni 


ES (^*(^.+i) 


Zijl3a,{sj)Ax'{xi{sj),yi{sj),Sj) 

V2 


x'iiSj + l) + 


V2 


Z£Ax'{xi{sj),yi{sj),Sj) 

V2 




( 20 ) 


where the coordinates Xi(sj+i),x'(sj+i) are those of 
the particle transported from Sj to Sj+i in presence 
of space charge, but without noise. In the following, 
we derive the artificial emittance growth averaged over 
many applications of space charge kicks. For this we 
hnd, that all terms proportional to Zij, or Zb have 
a minor contribution to the emittance growth, since 
they are in average zero with a variance proportional 
to the number of space charge kicks applied. The same 
argument holds true for the product of Zij and Zb, 
since they are statistically independent. The terms pro¬ 
portional to Zfj, and {Z'^j)'^ do not average to zero as 
{Zij)'^ = (Zb)^ = 1 V i,j, and thus they dominate the 
contribution to the emittance growth. All terms pro¬ 


portional to (Ax')^ are mnch smaller than any other 
term in Eq. 20, and can thus be neglected. Finally, as 
we consider matched beams, we use (x^) = (/l^x'^), and 
therefore we find 

ex(sj-i-i) — si{sj+i) + 

(x^(sj+i))((Ax'(x(sj),2/(sj),Sj))^). (21) 

Here, e^(sj+i) is the square of the emittance of the 
beam propagated without noise.The right hand side of 
this equation characterizes the averaged artificial emit¬ 
tance growth in one integration step; its explicit expres¬ 
sion is given by 
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A 4 (si+i) = {x^{sj+i)){i^x'{x{sj),y{sj),sj))^) = 


In the following, 5Ex{x, y, s) as given by Eq. 5 and Eq. 7, 
is used with the simplifying condition that the fluctua¬ 
tions due to the bi-linear interpolation are set to zero, 
i.e. ^[x,y^s) = 1. This ansatz well describes a fine 
grained mesh. 

The average of SE^(x, y, s) over x and y can be eval¬ 
uated by an integration over the particle distribution 
function, if the number of macro-particles Nm is large 
enough to approximate the p.d.f. Using the notation 
= ar{sj), we find for a Gaussian beam 


{5El{x,y,Sj)) = J dx J dyf{x,y)SEl{x,y,Sj) 


_ 1 

2 V Nm ) ■ 


Eor a K - V distribution we find respectively 

{dE^{x,y,Sj)) = - — -. 

Using Eq. 22, 23, and 24, we find in general 


Aelisj) = A 


V^/Vg f qSE^fi 


Nm ) x.j +1 


(23) 


(24) 


L(As)^ (25) 


with A a coefficient that incorporates the type of dis¬ 
tribution, with A = 1 for a KV, and A = 1/2 for 
a Gaussian. We thus find a dependence of the noise 
induced emittance growth on the number of macro¬ 
particles Nm, the number of grid-points Nq, the in¬ 
tegration length As and the variance tr^ of the dis¬ 
tribution. 

The factor of A = 1/2 for a Gaussian beam can be un¬ 
derstood in the following way: For a K - V particle dis¬ 
tribution, macro-particles are equally distributed, and 
the number of particles in each cell is therefore con¬ 
stant. In contrast to this, macro-particles are mostly 
located in the center for Gaussian beams. Therefore, 
the relative fluctuations of the electric field SE^/E^ are 
smaller, compared to a rms equivalent K - V beam. 

In this paper we intend to discuss the effect of PIG 
induced noise in a situation where the beam does not 
already exhibit emittance growth. Therefore, in the 
context of Eq. 21, we find that e^(sj+i), the squared 
emittance of the beam propagated without noise, will 
be equal to e^(sj ). In addition, the treatment of a con¬ 
stant focusing lattice further simplifies the formulas, as 
we find {x^{sj+i)) = = cr'^j- Therefore, if the 

beam has no intrinsic emittance growth, Eq. 21 becomes 

— slisj) + 

+{x^isj+i)){{Ax' ixisj),y{sj), Sj)f), (26) 


and Eq. 25 effectively provides the emittance growth 
Aex/A s due to PIG induced decorrelated noise: 


A£^ 

As 


f qSE,,o \\ 

- Am 2£AmocWj ^ ^ 


/ qAs 




2 

(x{sj+if){SEl{x{sj),y{sj),Sj)). 


( 22 ) 


B. Effect on the tune shift 


The artificial emittance growth may create another 
disadvantage in terms of an artificial change of tune- 
shift. In fact, space charge forces lead to a defocussing 
of the beam and thus to a change of the betatron tune 
Qx- For the K - V particle distribution (un-bunched) 
the associated tune shift is given by 


AQx{t) 


rpNp 1 1 

27r/?2y3 ’ 


(28) 


where tq = /{ATreomoc^) is the classical particle ra¬ 

dius. As the initial emittance Exit = 0) changes due to 
PIG induced noise, while the number of physical par¬ 
ticles Np remains constant, the associated tune shift 
changes by 


AQx{t) = AQxit = o f^^\.^\ (29) 

£x[t) 

By using the expression e^(t) = = 0) + Ae^, with 

slit) » Ae^, we find the approximation 


AQx{t) 


AQx{t = 0) 



Ael 

elit = oy 


(30) 


with the PIC induced emittance growth rate Ae^ as 
derived in the previous section. This derivation is valid 
for Gaussian beams as well, and the same scaling is 
found. However, the space charge induced tune shift is 
the double, see e.g. [27]. 


C. Application to the SISlOO 


In this section we benchmark the emittance growth 
predicted by Eq. 27 with simulations performed for the 
SISlOO heavy ion synchrotron [7] (without non-linear 
elements). We use the MICROMAP library [19] for 
tracking the particles. The PIC solver [16-18] makes 
use of a mesh of Nq x Nq = 64 x 64 grid points 
in a box of size Qax x day. For rms-emittances of 
£x = £y = 7.5 mm mrad the resulting space charge tune 
shift is AQx — —0.1686 and AQy ~ —0.1693, at tune 
position Qx — 18.87 and Qy ~ 18.72. The integration 
step As is set to 42 space charge kicks per betatron- 
wavelength. The results are presented in Fig. 4, where 
the red dots show the emittance growth from the sim¬ 
ulation, while the blue line is the theoretical prediction 
from Eq. 27, i.e. the random walk model applied to the 
full beam. 
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Number of turns 



FIG. 4: (Color) Evolution of the emittance e = 
due to numerical noise for a simulation with Nm = 10000 
macro-particles (top) and Nm = 20000 macro-particles 
(bottom) for a K - V beam tracked in the SISlOO lattice. 
The theory developed in this paper (blue line) predicts the 
artificial emittance growth in the simulation (red dots). The 
beam and simulation parameters are given in the text. 


For the simulations we use a coasting K - V beam. 
Therefore, we guarantee that space charge forces are 
linear and thus PIC induced numerical noise is the dom¬ 
inant cause of emittance growth. 

As seen in Fig. 4, the theory developed in this sec¬ 
tion recovers the artificial emittance growth found in the 
simulation of a complex machine. For decorrelated nu¬ 
merical noise the random walk model can thus be used 
in order to a) understand the effect of noise on single 
particles, and b) make predictions on the evolution of 
the whole ensemble of particles. 

Equation 27 can be used to find reasonable simulation 
parameters Nq, Nm and As. E.g. considering a sim¬ 
ulation with Nmi macro-particles, where an artificial 
emittance growth of Aex^Nmi i® observed, we predict 
that for a simulation with Nm 2 macro-particles 


A£x,Nm2 


Nmi 

Nm2 


A£x,Nmi ■ 


(31) 


The scaling law, Eq. 27, can also be used to find an 
optimal number of macro-particles for a different num¬ 
ber of grid points Nq. In fact, choosing the number of 


macro-particles such that 


= const, (32) 

Nm 


then the PIC induced emittance growth will be the same 
for all pairs No,Nm, hence one can choose the most 
efficient setting. 

Further, the agreement of theory and simulation al¬ 
lows to extrapolate the results for the SISlOO to a stor¬ 
age time of one second, corresponding to 2 • 10^ turns. 
We find, that for the simulation setup used above, we 
have to use Nm — 2.5 • 10® macro-particles, in order to 
limit artificial emittance growth to ~ 1%. If the physics 
case requires a higher resolution, such that we have to 
double the number of grid points Nq, we make use of 
Eq. 32. We then find, that to limit artifical emittance 
growth to ~ 1% for a storage time of one second, we 
have to use Nm — 3.5 • 10® macro-particles. A change 
of the integration length As changes the slope of artifi¬ 
cial emittance growth linearly. Thus, if we double As, 
we have to use the double amount of macro-particles 
to obtain the same amount of noise induced emittance 
growth. 

For the benchmarking, we considered space charge in¬ 
duced tune shifts smaller than the actual design goal for 
the SISlOO. We remark that for stronger space charge 
forces, the number of macro-particles has to be further 
increased, since the fluctuations become stronger, see 
Appendix C. A dedicated survey for the estimation of 
simulation parameters for certain operational scenarios 
will be part of future studies. 

Until now, only decorrelated numerical noise has been 
considered. However, the conjoint effect of the PIC al¬ 
gorithm and the particle tracker may cause correlation 
in the numerical noise. A study dedicated to this effect 
is presented in the following section. 


V. CORRELATED NUMERICAL NOISE 


In the previous section we discussed the effect of 
decorrelated numerical noise on single particles and the 
full beam. In this section, we study the case where nu¬ 
merical noise is correlated. If a coasting beam is tracked 
through a circular lattice, depending on the integration 
length As and the machine tunes Qx and Qy, the parti¬ 
cle ensemble returns (close) to its initial positions after 
a certain integer number m* of space charge kicks. Af¬ 
terwards, the same sequence of space charge kicks is ap¬ 
plied again, such as the electric field fluctuations. Thus, 
fluctuations of the electric field occur periodically and 
are no longer random. As the periodicity of the noise 
may create resonant effects, we call this a stochastic res¬ 
onance [28]. These kind of correlations may lead to an 
enhanced emittance growth that exceeds the estimates 
of the random walk model. 











A. Stochastic Resonances 

In mathematical terms, correlations in numerical 
noise are created if 

m* 

^A$j=27rn*, (33) 

where is the intra-kick phase advance in between 
Sj and Sj+i and n* is an integer. We define m* as the 
order of the stochastic resonance. 

For a constant focusing channel, where A$j = A<I> is 
constant for any j, Eq. 33 can be simplified to 

m*A$ = 27rn*, (34) 

with m*, and n* co-prime numbers. When this equation 
is nearly satisfied, the numerical noise on the first m* 
space charge kicks is completely random. However, for 
the following space charge kicks, the fluctuations of the 
electric field occur periodically according to 

6E^{xi, y*, Sj) ~ SE:c{xi,yi, Sj + m* As), (35) 

where, for simplicity, we assume the motion in y to be 
frozen. 

The recurrence of fluctuations has an important prop¬ 
erty if TO* is an even number: After to*/ 2 space charge 
kicks the phase advance is irn *, meaning that a particle 
i with initial coordinates Xi , yi is mirrored (antisymme¬ 
try) on the phase space axes to Xi,yi. In this case the 
periodicity relation Eq. 35 is modified, and reads 

5Ex{xi,yi,Sj) ~ -5Ea, (^i,y^,Sj + • (36) 

We observe that the recurrence of the electric field fluc¬ 
tuations described by Eq. 35 and Eq. 36 is limited to 
a maximum number of applications R of the m* space 
charge kicks. Afterwards correlations are canceled by 
accumulated noise on the particle’s positions. The per¬ 
sistence of correlations R depends on to* and on the 
granularity of the distribution that is controlled by Nq, 
Nm and tr,-. Therefore, a scaling law for artificial emit- 
tance growth in the presence of a stochastic resonance, 
analogous to Eq. 25 for decorrelated numerical noise, 
can only be derived by a dedicated study on the damp¬ 
ing of correlations, which is beyond the scope of this 
paper. 

The effect of the first to* kicks can be modeled by 
a random walk process, as described in Sec. IV A. The 
order of the stochastic resonance to* fixes the number of 
random elements in Zij and Z[y Thus, for a total of Nk 
space charge kicks in a simulation, if to* > Nk we re¬ 
trieve the random walk. If instead to* is a small natural 
number, correlations in the numerical noise are created. 
Depending on the number of kicks Nk, the persistence 
of the correlations R, and the order to* , correlated noise 
may enhance artificial emittance growth. This issue is 
discussed in the context of a periodic random walk in 
Appendix B. 

In the following, we outline the procedure for the two 
types of fluctuations described by Eq. 35, and Eq. 36. 


We start by considering numerical noise with a corre¬ 
lation as described by Eq. 35. We define two random 
variables Zij,Z[p for 1 < i < Nm, 1 < j < m*, in 
which Zij = ±1, and Zb = ±1 with equal probability. 
To model the effect of correlated numerical noise, the 
random elements Zij and Zb are applied periodically 
i?-times on all particles. 

When TO* is even, to model the noise with the prop¬ 
erty of Eq. 36, we define two random variables Zij, Zb, 
for f < i < Nm, 1 < j < m* 12, in which Zij = ± 1 , and 
Zb = ±1 with equal probability. In order to model the 
anti-symmetric property of the fluctuations, the next 
to */2 random numbers are defined by Zj = —'Z^ij 

and = —Zb, with f < j < vn* 12. This model 

is directly suggested by Eq. A5, observing that after 
m* 12 kicks the sign of aij, and bij inverts, see Ap¬ 
pendix A. 

Therefore, to* kicks of a resonance of an even or¬ 
der TO* have a similar effect on a particle as the first 
to*/ 2 kicks applied twice. This argument suggests that 
a resonance of order to */2 odd will produce a similar 
emittance growth as a resonance of order to* . 


B. Effect of correlations on the emittance growth 

In the following, we study the properties of the corre¬ 
lations in electric field fluctuations. For this purpose, a 
coasting beam is tracked in a constant focusing channel 
with machine tunes Qx — 0.19 and Qy ~ 0.24. The 
length of the channel is A = 3 m. For simplicity, we 
chose the integration length As = L. The electric field 
Ex{x,y,s) and Ey(x,y,s) at a: = y = 0 is calculated 
each turn for a total of 1024 turns, and a spectrum is 
retrieved by performing a Fourier transform of the data. 
In the corresponding spectrum, see Fig. 5, the highest 
peak occurs at the machine tune position. Other fre¬ 
quencies are excited due to a coupling of both planes. 

Since we pursue to study the effect of stochastic res¬ 
onances in one plane, we have to ensure that we do 
not have any interference of the other plane. This can 
be guaranteed by setting one of the tunes close to an 
integer value. If we choose Qx — 1.001, we decouple 
sufficiently the x and y planes, see Fig. 6 . 

In order to study the effect of the intra-kick phase ad¬ 
vance on the artificial emittance growth, we use again 
the constant focusing channel with the same simula¬ 
tion parameters as above. We now systematically scan 
the vertical tune Qy in the range Qy = 0.1...0.4, while 
keeping Qx — 1.001 fixed. The initial emittance in 
the y-plane is chosen for any machine tune such that 
we have a round beam {ax = cry). For each machine 
tune a coasting beam is tracked for 60000 turns and the 
slope of emittance growth per meter is calculated. The 
dependence of emittance growth on the machine tune, 
and thus on the excitation of stochastic resonances, is 
retrieved and shown in Fig. 7. 

By studying the simplified scenario of a constant fo¬ 
cusing channel, we guarantee a constant phase advance 
in between space charge kicks. Then, stochastic res¬ 
onances are excited when a multiple of the intra-kick 
phase advance is an integer of 2tt, see Eq. 34. For a re- 


9 



1 

0.9 

0.8 


g. 0.6 

c/5 

^ 0-5 

N 

^ 0.4 

E 

o 0.3 
z 

0.2 
0.1 
0 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Frequency 


FIG. 5: (Color) Spectrum analysis of the electric field in 
^-direction (top) and //-direction (bottom) aX x = y = Q. 
For this analysis, the beam is tracked in a constant focusing 
channel with tune Qx ~ 0.19 and Qy = 0.24. Peaks occur at 
the tune frequencies, while a coupling of both planes excites 
many other smaller local peaks. 



FIG. 6: Spectrum analysis of the electric field in //-direction 
at a; = // = 0. The vertical tune is Qy = 0.267 while the 
horizontal tune Qx — 1.001 is close to an integer value. A 
peak occurs at the tune position, while other frequencies are 
no longer excited. 
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alistic model of a machine the scenario gets more com¬ 
plex, as described by Eq. 33. The emittance growth 
rates, as presented in Fig. 7, confirm the prediction 
given by the theory developed in this paper. In the ab¬ 
sence of stochastic resonances (decorrelated numerical 
noise), we find a growth rate that can be predicted by 
the random walk model. The emittance growth rate has 
a (local) peak whenever a stochastic resonance is hit. 
The lower the order of the resonance m*, the stronger 
we see its effect on the emittance growth rate. Further, 
we find that if m*/2 is an odd integer, the growth rate 
is similar to the one of m*, e.g. in the case m* = 6. 

Until now, the space charge tuneshift was fixed for 
each simulation. In Appendix C we find that the width 
of the stopband of a stochastic resonance and the slope 
of artificial emittance growth are increased for more in¬ 
tense beams. 


C. Change of integration length 

For long-term tracking studies it is desirable to avoid 
the occurrence of correlated noise to a) minimize ar¬ 
tificial emittance growth and b) avoid uncontrollable 
artifacts, like tails in the distribution. In the following. 


FIG. 7: (Color) Emittance growth rate for different verti¬ 
cal tunes Qy, while Qx — 1.001. The tune Qy is changed 
systematically in steps of AQ = 0.001. 

we present a strategy to ensure that the numerical noise 
is decorrelated and thus controllable by a proper choice 
of simulation parameters, as discussed above. 

In most simulation studies the machine tune Qy is a 
multiple of the intra-kick phase advance, i.e. 

m 

= 2^gy, (37) 

i=l 

because the integration length As is a fraction of the 
machine length L. For simulations close to a machine 
resonance, i.e. of type nQy = A € N for n € N, stochas¬ 
tic resonances due to PIC noise will then be excited si¬ 
multaneously and will afflict the simulations. However, 
this situation can be avoided by choosing an appropri¬ 
ate integration length As. As a demonstration of the 
influence of the integration length on the occurrence of 
stochastic resonances, the study presented in Sec. VB 
is repeated for an integration length As = L ■ 0.95. The 
results are presented in Fig. 8. 

The blue line in Fig. 8 is the same as in Fig. 7, where 
As = L, while the red line corresponds to As = L-0.95. 
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Tune Qy 


FIG. 8: (Color) Emittance growth rate for different vertical 
tunes Qy, while Qx — 1.001. The tune Qy is changed sys¬ 
tematically in steps of AQ = 0.001.. The integration length 
is changed, compared to the blue reference line, by a factor 
of 0.95 (red line). 


strategy to avoid or minimize correlations in nnmerical 
noise, hence avoiding stochastic resonances. The theory 
developed in this paper can thus be applied to control 
unwanted noise effects in tracking simulations of high 
intensity beams. 

Last, we remark that in our studies we decoupled the 
PIC induced numerical noise in the transverse planes 
by setting the horizontal tune Qx close to an integer 
value. However, in more realistic applications, the ef¬ 
fect of coupling cannot be neglected, and certainly will 
introduce corrections. Investigations on the coupling 
between the planes will be part of future studies. 
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The peaks due to stochastic resonances are shifted for 
the modified integration length. 

The effect of changing the integration length can most 
easily be understood for the simplified case of a constant 
focusing channel. In fact, if the integration length is 
set to a fraction / of the initial one, i.e. As = L ■ f, 
the intra-kick phase advance is changed by a factor /. 
Therefore, the resonance condition, Eq. 34, is fulfilled 
for another tune Qy given by 

fl* 

m*Q'y = j . (38) 

We can thus effectively avoid correlated noise in our 
working point Qy by changing the integration length, 
and hence guarantee that single particles are solely af¬ 
fected by white noise. Then, artificial noise effects can 
be mitigated for an optimal parameter setting obtained 
using Eq. 27. 


VI. CONCLUSION 

In this paper we developed a model for predicting 
the effect of white noise on a particle beam by using an 
equivalent effective field fluctuation acting on all phase 
space coordinates. We use simulation results to infer 
the mathematical modeling of electric field fluctuations 
on a single particle in the beam, and use this model to 
develop a treatment of the full beam ensemble. The in¬ 
ferred scaling law as a function of the simulation param¬ 
eters enables us to estimate artificial emittance growth. 

We find that for specific machine tunes, which recall 
those of a machine resonance, PIC noise becomes cor¬ 
related. The effect of these correlations on emittance 
growth is more dramatic, and we explain this unusual 
dynamics in terms of stochastic resonances. Our studies 
for the constant focusing channel show that the intra¬ 
kick phase advance plays a crucial role on numerical 
noise as it defines the tunes at which stochastic reso¬ 
nances appear. This understanding is used to develop a 


Appendix A: Effect of the intra-kick phase advance 
on numerical noise 


In the following, the combined effect of linear track¬ 
ing and numerical noise is studied. For simplicity, we 
discuss the constant focusing channel. 

Let us consider a single particle initialized at position 
Xi,x^, and track it inclnding the space charge through 
a constant focusing channel. Due to numerical noise, a 
fluctuation Ax'{xi, x\, Sj) is induced at any longitudinal 
position Sj, where a space charge kick is applied. In 
total Nk space charge kicks are applied, i.e. the particles 
are tracked for a distance d = AsNk- The numerical 
noise at position Sj affects the particle’s phase space 
coordinates at sn^ by 




0 

ZUAx'{xi,x'i, 



(Al) 


We call the first component of this vector Oy and the 
second bij, which can be interpreted as random kicks 
on Xi resp. x'^ due to numerical noise at s = sj. Since 
we specialize on a constant focusing channel, we can de¬ 
compose the transfer matrix M into M = T~^R{A^)T 
with 


r= ( 

0 


and 


i?(A$) = 


cos(A<I>) sin(A$) \ 
— sin(A<I’) cos(A$) j 


(A2) 


(A3) 


Using I = Nk — j, we simplify Eq. Al to 

\Z'jAx'{xi,x'i,Sj) J \b. 

The random elements Oy, bij can be described by 
atj = Z[j sin (A$/) jdxAx'{xi,x'i, Sj), 
bij = Z[, cos (A^l) Ax'{xi,x'i,Sj). 


(A4) 


(A5) 
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A test-particle initialized at Xi{so), x'^{so) and tracked 
through the constant focusing channel accumulates the 
random elements , bij. For Nk space charge kicks the 
particle coordinates are 


Nk 

i,^Nk ) — ^^Nk ) A ^ ^ Ciij — Xii^S]^^ ) “k ? 

(A6) 

Nk ^ ' 

x'iisNk) = x'iisNk) + '^bij = Xi{sNk) + bi, 
i=i 

where Xi{sNk)yx'^isNk) ^^re the final coordinates of the 
test particle for a noise-free tracking. The random el¬ 
ements Qij, bij, as derived in Eq. A5, depend on the 
intra-kick phase advance A$. For most intra-kick phase 
advances A$ the cumulated numerical noise remains 
small, since the contributions a^- resp. bij are averaged 
to zero, i.e. (oi) ~ 0 and resp. (bi) « 0. We note that 
Qi (resp. bi) is a new random variable, which is the re¬ 
sult of the summation of the Nk random variables 
(resp. bij). The random variable ai has a variance of 


Nk 

2 _ 2 


(A7) 


because the random variables aij are decorrelated. A 
direct summation of this equation yields 


G 


2 

di 


Nk 


-[ 2 


72 


Px^x' 


(A8) 


which retrieves the property of the random walk. 
Therefore, we can substitute each aij with the random 
kick 


Zij—P,^Ax', 


(A9) 


with Zij = ±1 random, to retrieve the same final sta¬ 
tistical properties after Nk kicks. The same argument 
applied to the quantity bi shows that each bij can be 
substituted with 


bij —t Z[j-^Ax', (AlO) 

with Z'^ = ±1 random and statistically independent 
from the set Zij. This explains the effective random 
walk model used in Eq. 13. 

The random walk approach is not applicable any¬ 
more, if 


m*Ad> = 27rn*, (All) 

because a correlation is created. We show the breaking 
down of the random walk approach with a numerical 
experiment. For this, we define a sequence of random 
numbers Zij with the correlation 


Z,j = 7,(,+4) 1 < j < Afe - 4, (A12) 

which resembles correlated numerical noise in the pres¬ 
ence of a stochastic resonance. The cumulated effect for 
a constant fluctuation Ax' for Nk = 1000 is given by Oi 
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FIG. 9: (Color) Averaged cumulated numerical noise (oi) for 
100 initializations of Z'ij following Eq. A12. Peaks appear 
at phase advances A4> = 7r/2, and A'F = 37r/2. 
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FIG. 10: (Color) Averaged cumulated numerical noise (&i) 
for 100 initializations of ZP following Eq. A12. Peaks appear 
at phase advances A<1? = 0.0, A'F = 7r/2, A"!? = tt, A<1? = 
37r/2, and A4> = 27r. 







7t/2 


37t/2 271 


A<I> 



and bi, as described in Eq. A6. The averaged cumulated 
numerical noise (ui) and (bi) for 100 random configura¬ 
tions of Zb following Eq. A12, is shown in Fig. 9 and 
Fig. 10 for intra-kick phase advances 0 < A$ < 27r. 

We observe the cumulation of numerical noise for cer¬ 
tain intra-kick phase advances, i.e. |(ai)| >> 0, and 
|(&i)| >> 0. These results can be interpreted in the fol¬ 
lowing way. If the condition for a stochastic resonance, 
see Eq. All, is not fulfilled, then the elements a^-, bij 
are random and thus the random walk model can be ap¬ 
plied. Since both parameters are decoupled and are in 
average equivalent, the model given by Eq. 13 is justi¬ 
fied. However, in the presence of correlations this model 
is not valid anymore and the noise has to be treated dif¬ 
ferently, see Sec. V. 


Appendix B: Periodic random walk 

To explain the effect of a periodic random walk, we 
study the following simplihed case. A single particle 
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with the one dimensional coordinate x is affected by Nk 
random kicks of strength Ax = 1 in arbitrary units. We 
further enforce a periodicity of m*, i.e. we define a set 
of random numbers Xi G {—1,1} with i G {1, 2,.. .m*} 
that is applied H times. The number of repetitions i? is 
chosen to fulfill Nx = rn*R, with Nx constant. Since 
the averaged position after many random kicks is (x) ~ 
0 , the averaged second order moment can be associated 
with the variance. 

The dependence of (cr^) on m* is retrieved in a nu¬ 
merical experiment. For this purpose, is evaluated 
for 100 random initializations of Xi for each m* G 
{1, 2,... 2000}. The results are presented in Fig. 11. 



FIG. 11: (Color) Averaged second order moment (cr^) after 
Nk ~ 1000 random kicks for different periodicities m*. The 
blue line indicates the second order moment for the non¬ 
periodic random walk (without repetition). 

Fitting the numerical data, we find for m* < Nk the 
scaling of the averaged second order moment on m* 

(Bl) 

For m* > Nk, the averaged second order moment be¬ 
comes independent of m* since we apply only Nk kicks 
in this example. These results show that for a periodic 
random walk the diffusion is enhanced. The smaller the 
periodicity m*, the larger is the diffusion. 

When applying these results to the situation of a 
stochastic resonance created by PIC tracking, we have 
to consider that correlations in the numerical noise 
are washed out due to the numerical noise itself. The 
scaling, as given in Eq. Bl, may suggest that the 


periodic random walk could be an upper bound on 
the diffusion, because no loss of memory is considered. 
If the correlations are completely washed out, the 
diffusion type will the one of a random walk. This 
defines the lower bound on the diffusion, as indicated 
by the blue line in Fig 11. 

Appendix C: Effect of the strength of space charge 
on the excitation of stochastic resonances 

The excitation of stochastic resonances is certainly 
affected by the strength of the space charge force. To 
gain a quantitative knowledge, we repeat the study pre¬ 
sented in Sec. VB for different beam currents and thus 
different space charge forces. The artificial emittance 
growth is investigated in the vicinity of the third order 
stochastic resonance, which is presented in Fig. 12. 



Tune Qy 


FIG. 12: (Color) Emittance growth rate for different vertical 
tunes Qy in the vicinity of the third order stochastic reso¬ 
nance, while the horizontal tune is — 1.001. The vertical 
tune Qy is changed systematically in steps of AQ — 0.001. 
The colors represent different beam currents. 

As we learn from Fig. 12, an increase of the beam 
, and thus the space charge force, enhances artificial 
emittance growth and leads to a broadening of the stop- 
band of the stochastic resonance. The increase of ar¬ 
tificial emittance growth can be explained by the in¬ 
crease of the driving force of the resonance, i.e. of 
the space charge. At the same time, by increasing the 
space charge force we enlarge the tune spread of the 
beam. Therefore, the range of tunes Qy affected by the 
stochastic resonance grows larger. 
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